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NONLINEAR GRAVITATIONAL- WAVE MEMORY FROM BINARY BLACK HOLE MERGERS 
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ABSTRACT 

Some astrophysical sources of gravitational waves can produce a "memory effect," which causes a per- 
manent displacement of the test masses in a freely falling gravitational-wave detector. The Christodoulou 
memory is a particularly interesting nonlinear form of memory that arises from the gravitational-wave 
stress-energy tensor's contribution to the distant gravitational-wave field. This nonlinear memory con- 
tributes a nonoscillatory component to the gravitational-wave signal at leading (Newtonian-quadrupole) 
order in the waveform amplitude. Previous computations of the memory and its detectability considered 
only the inspiral phase of binary black hole coalescence. Using an "effective-one-body" (FOB) approach 
calibrated to numerical relativity simulations, as well as a simple fully analytic model, the Christodoulou 
memory is computed for the inspiral, merger, and ringdown. The memory will be very difficult to detect 
with ground-based interferometers, but is likely to be observable in supermassive black hole mergers with 
LISA out to redshifts z <2. Detection of the nonlinear memory could serve as an experimental test of the 
ability of gravity to "gravitate." 

Subject headings: black hole physics - gravitation - gravitational waves - relativity 



1. INTRODUCTION AND MOTIVATION 

In the typical picture of a coalescing binary black hole 
(BBH), the gravitational- wave (GW) signal oscillates about 
a zero value with an amplitude that grows slowly during 
the inspiral, rises rapidly during merger, and damps back to 
zero during the ringdown. This picture is not entirely accu- 
rate. In actuality, the GW signal displays a memory: for a 
circularized, inspiralling binary there is a growing, nonoscil- 
latory contribution to the"+" polarization that causes the 
signal amplitude to damp to a nonzero value (see the inset 
of Figure [1]). In a freely falling GW detector, this memory 
causes a permanent displacement of the test masses that 
persists after the GW has passed. While at late times the 
memory yields a constant (undetectable) shift in the space- 
time metric, the buildup of the memory is detectable. 

For gravitating systems with unbound components, a lin- 
ear memory effect has been known since the 1970s (see the 
references in Thorne 1992, hereafter Th92): this arises from 
nearly-zero- frequency changes in the time derivatives of the 
multipole moments of the source. For example, a hyperbolic 
binary has linear memory because the multipole moments' 
time derivatives depend on the (unbound) relative veloc- 
ity at late and early times, which has the same magnitude 
but changes direction. Linear memory also appears in sys- 
tems that undergo a kick (newborn neutron stars, recoil- 
ing black holes) or eject particles anisotropically (neutrino 
emission in supernovae, gamma-ray burst jets) [see Favata 
2009b (MFl) and Favata 2009a for references and further 
discussion] . 

I n the 1990's Blanchet fc Damour (1992, hereafter BD92) 
and lChristodouloul ()1991l ) independently discovered a non- 
linear memory effect that is present in all GW sources. This 
"Christodoulou memory" arises from the unbound gravi- 
tons radiated from the system (Th92): the lost GW energy 
contributes to the source's changing multipole moments. 
The Christodoulou memory is a unique manifestation of 
the nonlinearity of general relativity (GR) and is interest- 
ing for several reasons: (1) unlike other nonlinear effects 
on the GW amplitude such as "tails" (backscattering of 
GWs off of spacetime curvature), the nonlinear memory 
is nonoscillatory (it builds up over time as GW energy is 



lost); (2) even more so than tails, the memory is sensitive to 
the e ntire past history of the binary's motion (jArun et al.l 
120041 ): (3) even though the memory originates from the GW 
stress-energy tensor's contribution at 2.5 post-Newtonian 
(FN) and higher orders, the memory affects the GW ampli- 
tude at leading (Newtonian) order^ ; and (4) the memory is 
detectable and its observation could provide a strong-field 
test of GR. 

Previous calculations of the mem ory for circu larized bina- 
ries have treated onl y the inspiral (jWiseman fc Will 1 9911 : 
IBlanchet et al.l 120081 and have recently been extended to 
high PN orders (MFl). Estimates of the memory's de- 
tectability have either treated the memory as an unmodeled 
burst (Th92) or only included the buildup of memory dur- 
ing the inspiral (Kennefick 1994, hereafter Ken94). This 
work provides the first realistic estimate of the evolution 
of the memory, accounting for all phases of BBH coales- 
cence (inspiral, merger, ringdown). Properly modeling the 
memory during the merger and ringdown is especially im- 
portant because the reciprocal of the memory's rise time 
[- 1/(70M) ~ 1.45mHz(2 x WMq/M) for a BBH with 
total mass M and reduced-mass ratio rj = 0.25] can lie near 
the peak sensitivity of ground- and space-based GW de- 
tectors, affecting the power at those frequencies. Figures 
[1] and [2] illustrate the differences between previous memory 
computations and this work. 

Current numerical relativity (NR) simulations can best 
compute the dominant I = m = 2 mode of the waveform, as 
well as other (higher-order) oscillatory modes that contain 
no physical memory. Because these simulations cannot yet 
accurately extract the m = modes that contain memory, 
a semi-analytic calculation of this effect is necessary. (See 
MFl for further discussion.) 



^ Partly because of its high PN origin, partly because its nonoscil- 
latory nature makes its detection difficult, and partly because the 
effect is not widely known, the memory is usually ignored in standard 
discussions of the leading quadrupole-order waveform. 
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Fig. 1. — Evolution and saturation of tiic memory near tlic merger. 
The main plot shows several calculations of [Equation ^]. 

The solid (black) line uses the full EOB formalism calibrated to 
numerical relativity simulations (see the text). The dashed-dotted 
(green) line uses the same formalism but without any EOB amplitude 

corrections [Equation ((TOj with r^U -> {Mn)'^/'\ ^22/^2*^*^ = 1]. The 
dotted (cyan) curve is the minimal- waveform model [MWM; Equation 
^ with ^^max — 2 and rm — 3A4]. The short-dashed (blue) curve is 
the minimal-waveform model multiplied by a "fudge factor" 0.77 
so that it matches the fuU-EOB curve at late times. The long-dashed 
(red) curve is the inspiral memory truncated at rm — * 5M [Equation 
^ without the sum; this is the model used in Ken94] . All curves (ex- 
cept the "fudge factor" one) approach the same value at early times. 
The inset shows the h-^. waveform with (blue, solid) and without (red, 
dashed) memory computed using the fuU-EOB model [the oscillatory 
term s contain only the I = itm = 2 modes via Equations JlJ, ([Sj, and 
mOl ll. All plots are for equal-mass mergers with the matching to the 
ringdown signal near tm/M ^ 3522. 



2. CALCULATING THE MEMORY 

We begin by expanding the GW and "x" polariza- 
tions on a basis of spin- weighted spherical harmonics, 

+00 / 

V2R{h+-ih-^) = Y, J2 iUim~iVi„,)-2Y'"\e,<f). (1) 

1=2 m=-l 

Here (i?, 0, $) are the spherical coordinates pointing from 
the source to the observer; Uim and Vim are the coefficients 
of the radiative mass and current multipole moment ten- 
sors Z/^iii2...i, and Viiia-.-i; expanded on the basis of spherical 
harmonic tensors.^ The radiative moments are functions of 
retarded time Tr and appear in the general outgoing so- 
lution of the linear vacuum-wa ve-equation for GWs. The 
PN wave-generation formalism (jBlanchetl 120061 ) provides a 
nonlinear, iterative algorithm that relates these radiative 
moments to a family of source multipole moments that are 
expressible as integrals over the stress-energy pseudotensor 
of the matter and gravitational fields of the source. This 
iterative algorithm leads to the following weak-field expan- 
sion for the radiative mass multipoles (BD92; when not 
explicit, G = c = 1): 
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Fig. 2. — Illustration of how the memory model affects the signal- 
to- noise ratio (SNR). The short-dashed (blue) curve (the main result) 
shows the characteristic signal strength he of the total coalescence 
memory vs. frequency [computed us ing the minim al- w aveform model 
times the factor 0.77 and Equations I I12I I. JSjl, and I I14I I; see the associ- 
ated time-domain curve in Figure [l]. The SNR corresponding to this 
model is 8.9. Contrast this with the long-dashed (red) curve, which 
includes only the inspiral memory truncated at the last stable orbit 
(LSO) [see the associated time-domain curv e in Figure [T] the Fourier 
transform in this case is given by Equation I I14I I without the sum and 
with Tm — > 5Af]. The corresponding SNR, is 0.58. At low frequencies 
these curves approach the two horizontal (dotted) lines, which rep- 
resent a signal model that treats the memory as a step function in 
time, with the size of the step corresponding to the curves 

in Figure [1] that asymptote to 12.2 or 5.0. The corresponding SNRs 
are 23 and 9.5 for the upper and lower dotted horizonal lines. The 
thick solid (black) curve shows the LISA sky-averaged noise spectrum 
hn{f)- For comparison the he corresponding to the primary (oscil- 
latory) waves is also shown [the top solid (orange) line], truncated 
at the LSO frequency /lso = [5^^^-!rM{l + z)]-'^ . The SNR for the 
inspiral waves is 820. 



piece of the radiative multipole is given by (BD92, MFl) 



U, 



(mom) _ 32tT / {I ~2)\ 



Ira 



r.2-l 



2{l + 2)\ 



Tr 



dt 



dn^^ 

dtdVl 



(3) 

where ^^J^ is the GW energy flux. We ignore the radiative 
current moments, which have no nonlinear memory. 

For an adequate first estimate of the total coales- 
cence memory we focus only on the dominant I — 2 
multipoles. Substituting the leading-order energy flux, 
dE.„ - (327r)-iE™,™'4^4^'-2>^^™2r2™', intoEqua- 



dtdfl 



tion ([3]) and performing the angular integrals gives the time 
derivative of the leading-order contribution to the memory 
piece of the radiative mass moments: 
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where we have specialized to orbits in the xy plane (/2±i = 
0) and ignored moments that enter at higher PN orders. 
The UiQ for odd-Z vanish, and those with I > 6 enter at 

higher PN orders. The Ui^'^™'' for m 7^ yield oscillatory 
terms at 2.5PN and higher orders that do not contribute to 
the memory. 

Substituting Equations ([4]) into Equation ([T|) gives the 
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binaries and standard choices of the polarization tensors] : 
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To model the evolution of the source-quadrupole moment 
we follow the "effectivc-one-body" (EOB) approach (see 
Damour 2008 for references) calibrated to the results of 
NR simulations. The EOB framework attempts to extend 
the range of validity of the PN equations of motion to the 
nonadiabatic "plunge" region. It relies on a mapping of the 
PN 2-body Hamiltonian to the Hamiltonian of a point mass 
in a deformed Schwarzschild metric. 

It is instructive to first implement a simple and entirely 
analytic model for the coalescence that tries to qualitatively 
capture most of the important physics while minimizing 
complexity. This minimal-waveform model (MWM) con- 
sists of matching the leading-order inspiral moments to a 
sum of quasi-normal modes (QNMs). During the inspiral 
the g"^ derivative of l2±2 is 



Tinsp(g) _ r, 
^2±2 — ^ 



r;Mr'^(T2iw)«e 
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where = 93 = (M/r"^)^/^ is the orbital frequency, r — 
Tm(X ~ T/t„Y/'^ is the orbital separation, Lp is the OPN- 
order orbital phase, t„ = (5/256)(M/r/)(r„/M)'*, T = t- 
tm, and is the orbital separation at the "matching time" 
tm- For t > tm the quadrupole- moment derivatives are 
modeled as a sum of ringdown QNMs: 
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22ne 
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where a,. 



, + T;^^„, with QNM angular frequencie s 



^ eq uencie s 

i^imn and damping times Timn given bv lBerti et al.l (|2006D . 
These QNMs depend on the final mass Mf and the di- 
mensionless spin parameter a/ of the BH merger remnant, 
which are d etermined by NR s imula tions (I used the fits to 77 
in Table I of lBuonanno et aT] (|2007[ )). The coefficients A;,„„ 
are determined by matching Equations ©-(H]) a,tt — t,n for 
2 < (g,p-h2) < n„iax + 2.^ 

Substituting these relations into Equation ^ (and using 
dt = dr/f for the t < tm integral) yields 
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where H{T) is the Heaviside function. 

In addition to this simple, analytic model, I also imple- 
ment a more realistic EOB description of the moments that 
closely follows the approach of Damour et al. (2008b, here- 
after DNJ), in which EOB waveforms are calibrated to Jena 
and Caltech/Cornell NR waveforms: D uring the inspiral 
we so lve the EOB equations of motion [Damour fc NagaJ 
(|2007l) : Equations (7)-(ll)] for r, ip, and the canonical mo- 
menta Pr*: and p^. For the "radial potential" A{r) we 
use the (1,4)-Pade resummation of DNJ Equation (6) with 
as = 25. In the p^ equation we use DNJ Equations (8)-(ll) 
and Damour & Nagar (2008, hereafter DNCIT) Equations 



^2±2 



(17)-(18). For the parameters orr and Wpoie we use the val- 
ues in DNJ Table II for 05 = 25. For initial conditions we 
use ro = 15M, ip^ = 0, and the "post-post-circular" condi- 
tions of Damour et al. (2008a, hereafter DNAEI) Equations 
(l)-(3) for [pr'TPif)- We replace Equation ([7|) with 

2y^(T2z)^A/3-.(r^J7)3.-4gT2.^^^^yNQC^ 

(10) 

where = r^'^/'^, tp = [DNJ Equation (11)], n = (p, F22 

is given in DNCIT Equations (5)-(ll), and f^^^ = [DNJ 
Equation (12)] (with 6 = and a given by the linear fit 
in DNJ Sec. III). For the ringdown we use Equation ([H]) 
above with five QNMs (with modes chosen as in DNAEI, 
Section III). To determine the coefficients A22n we match 
f 2) . 

I22 at five points centered around the time when the ra- 
dius equals the EOB light-ring [the peak of ^{t)] as de- 
scribed in DNJ (to improve the fit I shift the matching 

time by — 3Af). Once constructed 122'' substituted into 
Equation ([6]) and numerically integrated using the initial 

value /i("""'(0) = SirM/ro. The results of the EOB and 
the minimal- waveform models are shown in Figure [1] Note 
that the EOB corrections reduce the memory's magnitude. 

3. SIGNAL-TO-NOISE RATIOS FOR THE MEMORY 

To compute the memory's signal-to-noise ratio (SNR) we 
use the MWM [Equation dH)] multiplied by a "fudge factor" 
~ 0.77 that ensures late-time agreement with the fuU-EOB 
model (see Figured]). The sky-averaged SNR^ for a detector 
with single-sided noise spectral density S'„(/) is 



(SNR2 



hUf) df 
hlif) f ■ 



(11) 



Here hn{f) — yaJS^Jj) is the sky-averaged rms noise 
amplitude, where a ~ b for o rthogonal a rm detectors like 
LIGO and a = 20/3 for LISA (jBarack fc C utler 200^). The 
memory's characteristic amplitude is defined by 



Kif) = 2{l+z)f{\hf''-^>[{l + z)f]\')y-'\j,^oU(i+^)^ (12) 

where h!"^°"^\f) denotes the Fourier transform (FT) of 
Equation ([5]), Dl{z) is the luminosity distance, and we use 
the cosmological parameters ifo/(100 km/s/Mpc) « 0.70, 
= 0, ~ 0.28, and w 0.72. 
To compute the FT we use: (1) the FT of H{T), 
T[H{±T)] = [S{f)±i/{Trf)]/2; (2) t he FT of g(r)e-^^ for 
> ; and (3) Equation (13.2.6) of Abramowitz fc StegunI 
(|1972D for the FT of the first term of Equation ®: 



(1 - r/r„)i/4 



dr = r„f/(l,7/4,27ri/T„), (13) 



where U is Kummer's confiuent hypergeometric function of 
the second kind. For / > the FT of Equation ^ is then 



r (mem) 
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iTTJ I Tm 



— y 



' 22n" 22n' ^22n^22n' 



^^„tto2''*^^(^22"+^22„'). 



(14) 



Combining the above formulae with sensitivity curves for 
the various GW interferometers allows us to estimate the 
memory's detectability. Unless the source is within the 
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Fig. 3. — Angle- averaged signal-to-noise ratio of the memory for 
equal-mass LISA binaries as a function of total binary (source-frame) 
mass M and redshift z. A low-frequency cutoff of 10^^ Hz in Equation 
JTTJl was assumed. 

Local Group, the memory signal will be too weak to be 
detected with current LIGO. Advanced LIGO will have 
a ten times greater sensitivity, yielding a SNR « 8 for a 
50Mq/50Mo BBH at 20 Mpc. For LISA the prospects of 
detecting the memory from supermassive BBH mergers are 
much better: a IO^Mq/IO^Mq merger at z = 2 has a SNR 
for the memory of « 8.9 (see Figure^]). Figure [3] shows the 
SNR as a function of mass for equal-mass LISA binaries at 
selected redshifts. 

Figure [5] also illustrates the sensitivity of the SNR to the 
memory model. Ignoring the merger and ringdown signif- 
icantly underestimates the memory's SNR. Modeling the 
memory as a step function in time over-estimates the mem- 
ory's SNR, even if the correct saturation value A/i*^™"™^ is 
used. This approximation is implicit in the SNR estimate 
of Th92 and is equivalent to using the zero-frequency-limit 
f^i"r'(.f ) ~ iA/t('^''" ')/(27r/)] to approximate the mem- 
ory's FT (|Smardll977n . 

4. DISCUSSION AND CONCLUSIONS 

The Christodoulou memory is a dramatic example of how 
the nonlinearities of general relativity can manifest them- 
selves in a detectable GW signal. This first attempt at 
including the merger and ringdown illustrates the impor- 
tance of the full coalescence in any model of the memory. 



This calculation also indicates that despite the recent suc- 
cesses of numerical relativity (NR), analytic methods can 
still be useful in understanding the highly nonlinear regime 
of BH mergers. It is especially interesting to see that a sim- 
ple analytic description like the "minimal-waveform model" 
can qualitatively describe the memory. This is partly due 
to the memory's independence of the GW phase, which is 
more sensitive to PN corrections than the amplitude [see 
Equation ©]. 

Although based on an FOB model calibrated to NR sim- 
ulations, these memory estimates relied on various approx- 
imations. While the memory is not easily extracted from 
current NR simulations (MFl), input from NR is needed 
to compute the memory accurately. Preliminary results 
from a hybrid PN/NR calculatio n using the I = m = 2 
Caltech/Cornell merger waveform (jScheel et al.ll2009l ) indi- 
cate that the memory saturates at a value A/i('°°™^ « 9.6 
(M. Favata 2009, in preparation). This « 27% difference 
with the fuU-EOB model used here [A/i^'^'^™) « 12.2; see 
Figure [1] is roughly consistent with the « 20% error be- 
tween the amplitude of the NR waveform and the FOB 
model of DNJ near merger (see their Figures 8 and 9). Fu- 
ture work will consider higher multipole in teractions and 
recen t improvements in the FOB formalism (jPamour et al.l 
l2009f ). 

Signal-to-noise ratio (SNR) calculations suggest that the 
memory from stellar-mass BH mergers is unlikely to be de- 
tected with Advanced LIGO, but the memory from super- 
massive BH mergers is potentially detectable by LISA out 
to a redshift z < 2 (Figure[21). This detectability assessment 
is based on the memory's SNR lying above some threshold 
~ 5-8. A better measure of detectability should consider 
separating the memory from the stronger oscillatory merger 
and ringdown waves. This will be add ressed in later work . 
Future Mock LISA Data Ghallenges (jArnaud et al.ll2007f) 
should consider incorporating the memory into their data 
sets. The memory also introduces new parameter depen- 
dencies into the GW signal, and it would be interesting to 
see how th is affects the estimation of binary parameters 
()Arun et a l. 2008). If the nonlinear memory is eventually 
detected, it could provide an interesting test of general rel- 
ativity, particularly of the ability of gravity to "gravitate." 
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